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Artificial magnetic fields (AMFs) created for ultra cold systems depend sensitively on the internal 
structure of the atoms. In a mixture, each component experiences a different AMF depending on 
its internal state. This enables the study of Bardeen-Cooper-Schrieffer pairing of fermions with 
unequal effective charges. In this Letter, we investigate the superconducting (SC) transition of 
a system formed by such pairs as a function of field strength. We consider a homogeneous two- 
component Fermi gas of unequal effective charges but equal densities with attractive interactions. 

We find that the phase diagram is altered drastically compared to the usual balanced charge case. 

First, for some AMFs there is no SC transition and isolated SC phases are formed, reflecting the 
discrete Landau level (LL) structure. SC phases become reentrant both in AMF and temperature. 

For extremely high fields where both components are confined to their lowest LLs, the effect of the 
charge imbalance is suppressed. Charge asymmetry reduces the critical temperature even in the 
low-field semiclassical regime. We discuss a pair breaking mechanism due to the unequal Lorentz 
forces acting on the components of the Cooper pairs to identify the underlying physics. 

PACS numbers: 


Cold atom experiments have realized novel many- 
particle systems, challenging some of the most funda¬ 
mental models of condensed matter theory. In particular, 
Bardeen-Cooper-Schrieffer (BCS) theory of fermion pair¬ 
ing, which has successfully explained superconducting 
(SC) and superfluid (SF) behavior in a large number of 
systems, had to be extended to cover new regimes. Pair¬ 
ing due to resonant interactions has been explored both 
theoretically and experimentally, uncovering the BEC- 
BCS crossover in detail. Density imbalance between com¬ 
ponents forming the Cooper pairs and resulting uncon¬ 
ventional SC states were first considered for condensed 
matter systems, but, experimental observation of polar¬ 
ized SFs [H, i with cold atoms required significant im¬ 
provements upon prior approaches. Similarly, Cooper 
pairs made up of fermions with unequal masses have been 
explored theoretically ii- 

The constituents of cold atom experiments are neutral 
atoms. The dominant interaction between these atoms 
are through s-wave scattering which can be tuned via 
Feshbach resonances between the atoms. While the ab¬ 
sence of Coulomb interactions facilitated the realization 
of some fundamental condensed matter models, the neu¬ 
trality of the particles prevented the observation of the 
effects of an external magnetic field on these systems. 
Initial efforts in this direction used rotation to mimic the 
magnetic field which brings further constraints on the 
confining potential of the ultracold system 

Over the last five years, a significant development, 
namely the creation of Raman laser-assisted artificial 
magnetic fields (AMFs) for neutral atoms [ 1 ], has ex¬ 
tended the capabilities of cold atom experiments. These 
AMFs are realized by coupling the internal states of the 
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atoms to light to imprint a Berry phase on the motion. 
While a number of different schemes have been used to 
manufacture these synthetic Hamiltonians, all of them 
sensitively depend on the internal excitation structure 
0 . Hence, for a mixture of two different atom species or 
even a mixture composed of atoms in different hyperfine 
states, the effective magnetic field acting on each com¬ 
ponent can be different. For example, the ^-factors for 
®'^Rb 5 S'i/2 a = 1 and 5S'i/2 F = 2 have a 3/2 
ratio. If the scheme in Ref.@ is applied to a mixture of 
these atoms, position dependent detunings, consequently 
the AMFs, would reflect this ratio. Although the Zeeman 
shifts due to the real magnetic field are utilized to create 
the AMF, this artificial field only couples to the spatial 
motion of the atoms and does not cause an artificial Zee- 
man effect Q- 

In this Letter, we explore the consequences of an AMF 
that couples unequally to the fermions forming a Cooper 
pair. Essentially, we consider the pairing of fermions with 
different cyclotron frequencies, which we regard as un¬ 
equal effective charges coupling to the same AMF. We 
discuss the conditions for pairing and the response of the 
paired state to the external AMF. We show that this sys¬ 
tem displays reentrant SC in temperature, i.e. a normal 
sample at zero temperature can become SC as temper¬ 
ature is increased. Oscillatory dependence of Tq on the 
AMF, which is a direct consequence of the Landau Level 
(LL) structure of single-particle excitations, is observed. 
However, for some AMFs, SC state is not preferred even 
at zero temperature. We calculate the phase diagram of 
the system for various representative cyclotron frequency 
ratios and present physical mechanisms to elucidate the 
fundamental changes in the SC transition. 

We consider a mixture of two fermion species of equal 
mass and equal density. The system is assumed to be 
spatially homogeneous, as in most cold atom experiments 
the effects of the confining potential can be taken into ac- 



2 


count through the local density approximation. An AMF 
of arbitrary strength is acting on the system by coupling 
only to the orbital motion of the fermions but causing no 
Zeeman shift. The coupling of the AMF to each compo¬ 
nent is different, defining the effective charges qi and 52 - 
Corresponding cyclotron frequencies wi = qiB/m and 
UJ 2 = q 2 B/m define the respective LL separations. We 
introduce the relative frequency cOr = a; 2 /wi and the ef¬ 
fective frequency w = Within the Landau gauge 

A = (0,i3a:, 0), the non-interacting Hamiltonian can be 
written as 


'Ho(wi,a;2) = '^£iufLfi^ + e2^fl^f2u: (1) 

V 


where the index v = (n,ky,kz) incorporates the LL in¬ 
dex n, momentum along the ^-direction kz, and the mo¬ 
mentum ky which also labels the LL degeneracy. The 
associated kinetic energy is = tiy — = h^k^/2m + 

— /ij, with i = 1,2. The chemical poten¬ 
tials /ii are not equal but are chosen to fix the den¬ 
sity of both species to be the same at each AMF value, 
Ni = N 2 = N/2 = Bqt/2TTh'^^ J dkzF{eiy) with the 
Fermi-Dirac distribution function F{eiy). The parti¬ 
cle densities ar e then scaled with the effective magnetic 
length, i = •Jhfmjjj, i.e. ni = We numeri¬ 

cally solve the number equations for chemical potentials 
at each AMF. Thus, for a fixed charge ratio uir and total 
real-space density N, changing the AMF strength alters 
only the effective density ni. 

This single particle spectrum is unique as the LLs of up 
and down spins do not match in energy. Since their zero 
point energy and the separation between the LLs are dif¬ 
ferent, two LLs can have equal energy only if the charge 
ratio ujr is a rational number. When the chemical po¬ 
tentials are adjusted to equate the densities, low energy 
single-particle excitation spectra for up and down spins 
are asymmetric. This mismatch has drastic consequences 
on pairing when the interactions are introduced. 

The two species interact resonantly through s-wave 
scattering which we model using the two-channel Hamil¬ 
tonian following Refs.@, @ studying the balanced mag¬ 
netic field case, uj^ = 1. We write the interacting Hamil¬ 
tonian, 


n{uJi,UJ2) = + £2yfLf2y + 

V 

+ aY,QiyfiJ2y,FH.c. ( 2 ) 

W! 


The open channel fermions interact to form the closed 
channel boson with energy £ 3 = 7 -!- (hioi + ^ 102 ){2 — 
k-i ~ ^-2 + C, where 7 is the unrenormalized detuning 
between the closed channel boson and the open channel 
fermions. C is the counterterm which is set to compen¬ 
sate the divergence of the boson self energy due to the 


infinite number of LLs of the closed channel fermions, 

(n-I-n')! 

AirfF nln'l (£? + 

X , ^ (3) 

Y TZ -|- ^ -|- UJr (n' + ^ ) 

The closed channel boson eigenstates have the same form 
as the fermion single particle states, however the ef¬ 
fective charge of the boson is qs = < 71 + 92 , and its 
mass 2m. C onseque ntly the boson magnetic length is 
£b = £\£ 2 l\/£\ + £2 = \/^/m(wi + a; 2 ). Two fermions 
with up and down spin in different LLs can interact to 
form a boson through the coupling constant a and the 
overlap integral Qyyi which has the form 


Quy' — ^{kz + kz') 




V^V^2"+"'n!n'!(£2 + £2)„+n'-si/2 

e-’^-^^^l+^^^Hn+u'(ky^£l+£tj. (4) 


where i?„(.) is the n-th Hermite polynomial. The dom¬ 
inant pairing mechanism is through the closed channel 
boson with the lowest possible energy. Hence, the wave 
function of the boson is in general a superposition over 
all the states in the bosonic lowest LL (LLL) with zero 
kz- The distribution of the boson wave function over the 
degenerate LLL states with different ky does not affect 
thyihysical properties after the bosons are integrated out 
Thus, we calculate the overlap integral only for the 
LLL state with ky = 0. The renormalization parameters 
7 and a are chosen in such a way that the system pro¬ 
duces low energy scattering properties for (a;i,a; 2 ) — t 0 . 
Hence, they are related to the physical parameters Og- 
scattering length and ro-effective range as follows. 


7 _ m 1 _ rpm^ 

Us/i^dTr’ 

In order to analyze the pairing of unequal charges, we 
examine this Hamiltonian around the SC transition. Fol¬ 
lowing mean-field theory, we introduce a < b >= A 
which is the average amplitude of the bosonic wave func¬ 
tion defining the order parameter. Near the transition, 
A is small and we expand the free energy. 


Pairing is favorable if the coefficient of the second order 
term is negative. The critical temperature for pairing is 
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FIG. 1. (Color online) Phase diagram of the system as a func¬ 
tion of dimensionless temperature T = ksT /hw and effective 
densit y m = where N is the real-space density and 

£ = Notice that increasing AMF, B = moj/y/qiq 2 , 

corresponds to lower m values. Two frequency ratios (Ur = 1 
and oJr = 0.95 are displayed for the same interactions strength 
ds = = 0.53. For Ur = 1, there is SC transition at 

any field, the oscillations in Tc (stars) originates from the LL 
structure. For ojr = 0.95, these oscillations evolve into bubble 
SC regions (shaded areas). The system is not SC even at zero 
temperature between the bubbles and the transition becomes 
weakly reentrant in temperature. At the low (many LLs) and 
high field (only LLL) regimes, Tc is not affected significantly 
by a small charge imbalance. 


obtained by setting this coefficient to zero 

_ J_ _h^^ (n-bn')! 

n- m nln'l (£?-I-£2)”+"’'+^ 

dfcz [ tanh ^ + tanh ^ 2 | 

27r ^ Ein + S2n' ^In + E2n' j 

(7) 

The chemical potentials are chosen so that the real- 
space densities of the two components are the same. The 
right-hand side of Eql3 the pairing susceptibility, deter¬ 
mines the behavior of Tq and can be used to understand 
the underlying physical picture for its evolution. 

We solve EqjT] numerically by calculating the pairing 
susceptibility for a given value of T and rii. The phase 
diagram of the system is then obtained by comparing this 
value with — l/og. For the numerical solution we scale all 
energies by the effective magnetic energy huj. Similarly 
the dimensionless scattering length is ds = 

Our equations are symmetric for Wr —^ which is equiv¬ 
alent to switching the indices of the components. We 
checked this symmetry numerically and concentrate on 
0 < Wr < 1 in the following. 

In FigITJ we present the phase diagram of w,. = 1 which 
is in agreement with Ref.Q. For the balanced charge 
case, there is always a critical temperature below which 
the SC state is preferred within the mean-field approxi¬ 
mation. The critical temperature is non-monotonic with 



the applied field. It first decreases and becomes expo¬ 
nentially small as smaller number of LLs are involved in 
the pairing, but, then increases when only the LLL con¬ 
tributes. This l^h-field SC has been studied for both 
solid state [13, [HI and cold atom systems Q . The os¬ 
cillatory behavior of Tc with the applied field is a direct 
result of the underlying LL spectrum. When the LLs of 
the two components coincide in energy, the pairing sus¬ 
ceptibility diverg es as IjVT at the peaks and as ln(T) 
at the rest in Fig 2(a)[ always guaranteeing a SC state. 


We also display the phase diagram for ujr = 0.95 in 
FigUl First of all, unlike the balanced case, there are 
AMF values for which a SC state is never favored. Oscil¬ 
latory behavior in pairing due to the LL structure causes 
stable islands of SC in the phase diagram isolated from 
the zero field SC phase. Although the low field SC is not 
affected from unequal frequencies as can be expected, the 
high-field SC proves to be surprisingly resilient too. The 
destruction of the low temperature SC in the presence 
of a small misalignment between LLs has been predicted 
[H, however, this misalignment also leads to a third ef¬ 
fect, reentrant SC. For some AMFs, at low temperatures 
the sample is in the normal state while at higher tem¬ 
peratures it becomes SC. Consequently, even for a slight 
asymmetry between the charges of the components, the 
phase diagram undergoes fundamental changes. These 
changes are most pronounced in the regime where only 
a few LLs are populated for both components and we 
display representative phase diagrams in Figl2] In the 
following, we discuss the physical reason for each of these 
three features. 


Fig 2(c) displays the pairing susceptibility for ojr = 
0.95 by focusing on intermediate field strength. The most 
striking feature of this phase diagram is the emergence of 
isolated islands of SC Fig |2(d)| which come about because 
Eql7]does not have a solution at some AMF strengths. 
The absence of solutions even for a minute amount of 
charge imbalance is best understood by considering the 
single-particle spectra. The one-particle density of states 
(DOS) for each component has sharp peaks at each LL 
threshold. For the balanced case, the DOS and the chem¬ 
ical potentials of the components are always equal. If the 
temperature is low enough, only the DOS near the chem¬ 
ical potential is relevant. At low temperatures, the pair¬ 
ing susceptibility diverges as at each LL threshold 

and the peaks in Tc follow from the one-particle DOS. 
When uJr y 1, LLs of different components do not have 
the same energy or total degeneracy. In general, chemi¬ 
cal potentials of both components must be chosen differ¬ 
ently to give equal real-space densities. If the mismatch 
between the chemical potentials and the LLs is large, the 
energy cost of exciting particles may not be redeemed by 
attractive interactions even at zero temperature. 


The pockets of SC phases roughly correspond to popu¬ 
lation of a new LL. Whenever a chemical potential crosses 
a LL threshold, there is a new set of states at the Fermi 
surface which become suddenly available to contribute to 
the pairing. The most favorable case for pairing is when 
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FIG. 2. (Color online) Pairing susceptibility (left panel) and 
respective phase diagrams (right panel) for three frequency ra¬ 
tios a;r = 1, 0.95, 0.75. The system is made dimensionless with 
effective magnetic energy tiw as in Fig[T] and decreasing ni 
corresponds to increasing AMF at fixed real-space density N. 
The phase diagrams are in linear scale in order to cover Tc = 0 
and plotted for intermediate held strengths where charge im¬ 
balance effects are most prominent, a) Pairing susceptibility 
of equal charges diverges at low temperature, guaranteeing 
SC for any held value. Divergence is more pronounced at 
LL thresholds. The corresponding phase diagram (b) is ob¬ 
tained by Eq[71 The phase boundary is also highlighted on 
the surface, c) Even a slight asymmetry between the charges, 
uJr = 0.95, lifts the low temperature divergences and the os¬ 
cillations in (b) turn into bubble SC phases (d). Each LL 
susceptibility peak is split into smaller peaks, thus, the bub¬ 
ble phases branch into smaller bubbles for weaker interac¬ 
tions (not displayed). e,f) The mismatch between LL spectra 
is greater for smaller uir, resulting in prominent reentrance 
with temperature. The maximum reentrance temperature is 
controlled by h\uj\ — uj 2 \- 


both chemical potentials simultaneously cross a new LL. 
For small charge imbalance, these threshold crossings 
happen within a small difference in AMF which essen¬ 
tially turns the Tc oscillations of the balanced case into 
bubble SC phases in Fig |2(d)| For a general Wr ratio, the 
picture is much more complicated. For chemical poten¬ 
tials to give equal densities and cross LL thresholds si¬ 
multaneously, UJr must be close to a simple fraction. This 
complicated behavior is evident in the pairing susceptibil¬ 
ities displayed in Figl^] where simple peaks of balanced 
LLs are split into smaller structures. For very strong 
attractive interactions, these bubble phases are not re¬ 


solved as Tc becomes comparable to the LL separation. 
A similar effect was predicted for the balanced case Q. 

As the AMF is increased further, we observe that 
Tc increases and reaches the same value with the equal 
charge limit. This surprising revival happens only when 
both components populate their LLLs. Although the de¬ 
generacies of the two LLLs are not equal, they both in¬ 
crease with increasing AMF. The chemical potentials of 
both components then lie very close to the correspond¬ 
ing LLL thresholds. Hence, the excitation cost for pairing 
decreases at such high fields. We can estimate the AMF 
for which the effect of charge imbalance vanishes by re¬ 
quiring all the particles with smaller effective charge to 
reside in their LLL. This estimate is in good agreement 
with our numerical results. 


Another fundamental change brought about by the 
charge imbalance is SC that is reentrant with temper¬ 
ature. While this effect is not clear for small imbalance 
as in Fig |2(d)| we found that it is a common feature of 
the phase diagrams for general uJr- A more prominent 
reentrant SC phase can be observed for uJr = 0.75 as 
in Fig 2(f) For some field strengths, the system prefers 
normal phase at zero temperature and becomes SC only 
above Tci- SC phase subsequently disappears after a 
higher temperature Tc 2 - Similar reentrant behavior was 
predicted for graphene bilayers [T^ and asymmetric nu¬ 
clear matter |14| . In our system it is easy to understand 
the physical basis for this reentrance. Increasing tem¬ 
perature generally prefers a disordered state, however, it 
also excites a significant amount of particles to higher 
LLs. Because of the antisymmetric and oscillatory na¬ 
ture of the DOS in the charge imbalanced system, states 
which are close to not only the chemical potential but also 
to LL thresholds are most favorable for pairing. Thus, 
if the pairing contribution from thermally excited parti¬ 
cles overcomes the entropy cost, increasing the temper¬ 
ature can drive the SC transition. With this scenario 
we expect the maximum lower critical temperature Tci 
to be of the order of LL mismatch between the compo¬ 
nents which agrees with our numerical results. In con¬ 
trast to other reentrant SC phases where a competing 
order precludes SC at low temperatures [l^ , the current 
system has reentrance solely due to non-trivial nature of 
the single-particle DOS. 

Although we concentrate on the interplay between dis¬ 
crete LL structure and charge asymmetry, it is worth 
mentioning that there is a profound effect even in the 
semiclassical regime where many LLs are filled for both 
components. The effect of an external magnetic field on 
a Cooper pair is modeled by only considering the phase 
acquired by the center of mass motion. This semiclassi¬ 
cal approximation due to orbital dephasing describes the 
upper critical field Hc 2 in type-II SC successfully. How¬ 
ever, if the charges of the fermions forming the pairs are 
different, the magnetic field couples the center of mass 
motion with the relative coordinate. As pairing is con¬ 
trolled by the relative coordinate, especially for tightly 
bound pairs, it is possible for the Lorentz force due to 
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the center of mass motion to break the pairs. Classically, 
if the center of mass of a bound pair with charges qi , q 2 
is moving with velocity u in a perpendicular magnetic 
field B, Lorentz force difference between the two parti¬ 
cles is F ~ Iqi — q 2 \Bv. This real space picture can be 
utilized to give a rough estimate for the strength of this 
pair breaking mechanism by comparing the work done by 
this force over the size of the pair to the SC gap. This 
effect becomes dominant especially for large charge ra¬ 
tios and we estimate the upper critical field due to this 
pair breaking mechanism as Hcs ~ Hc 2 which is 
in agreement with our numerical results. While compet¬ 
ing orders such as charge density waves, which were not 
taken into account in our approach, may complicate the 
physics in the high magnetic field limit, the decease of 
Tc with charge imbalance is observed even when many 
LLs are filled and the mean-field approximation is most 


reliable. 

In summary, the cold atom experiments with AMFs 
can create mixtures where each component has a different 
effective charge. The pairing between fermions of unequal 
effective charges presents a unique extension of BCS the¬ 
ory which is fundamental in diverse areas of physics. In 
this Letter, we find that even for a slight asymmetry be¬ 
tween the charges, the phase diagram changes drastically 
with the emergence of reentrant SC both in temperature 
and AMF. The oscillatory behavior of Tq with AMF 
for the balanced case modifies into isolated SC phases. 
For extremely high AMFs where both components are 
in their LLLs, the transition temperature is independent 
of the charge ratio. Finally, we argue that Tc is reduced 
due to pair breaking facilitated by unequal Lorentz forces 
on the charges forming the pairs. 
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